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ABSTRACT 

We use low-dispersion spectra obtained at the Magellan Observatory to study the broad Ha emission 
from the reverse shock of the infant supernova remnant SNR1987A. These spectra demonstrate that 
the spatio-kinematic structure of the reverse shock can be distinguished from that of the circumstellar 
ring and hotspots, even at ground-based spatial resolution. We measure a total dereddened Ha flux 
of 1.99(±0.22)xl0 -13 ergs s -1 cm~ 2 at an epoch 18.00 years after outburst. At 50 kpc, the total 
reverse shock luminosity in Ha is roughly 15 Lq, which implies a total flux of neutral hydrogen atoms 
across the reverse shock of 8.9xl0 46 s _1 , or roughly 2.3xl0 -3 Mq yr _1 . This represents an increase 
by a factor ~ 4 since 1997. Lyman continuum radiation from gas shocked by the forward blast wave 
can ionize neutral hydrogen atoms in the supernova debris before they reach the reverse shock. If the 
inward flux of ionizing photons exceeds the flux of hydrogen atoms approaching the reverse shock, 
this pre-ionization will shut off the broad Lya and Ha emission. The observed X-ray emission of 
SNR1987A implies that the ratio of ionizing flux to hydrogen atom flux across the reverse shock is 
presently at least 0.04. The X-ray emission is increasing much faster than the flux of atoms, and if 
these trends continue, we estimate that the broad Lya and Ha emission will vanish in <7 years. 
Subject headings: circumstellar matter — shock waves — supernovae: individual (SN1987A) — su- 
pernova remnants 



1. INTRODUCTION 

The collision between the ejecta of SN1987A and its 
circumstellar ring is now in full bloom, signaling the birth 
of the supernova remnant SNR1987A. This interaction 
was predicted (Luo & McCray 1991; Luo et al. 1994; 
Chevalier & Dwarkadas 1995; Borkowski et al. 1997) 
shortly after the discovery of the circumstellar ring. It 
began about a decade after the outburst with the dis- 
covery of the first of many "hot spots" (Sonneborn et al. 
1998; Michael et al. 2000; Pun et al. 2002). These hot 
spots are thought to result when the forward blast wave 
encounters and transmits radiative shocks into dense pro- 
trusions or "fingers" pointing to the interior of the ring. 
Since then, hot spots have encircled the entire ring (Sug- 
erman et al. 2002). 

Behind the blast wave, the expanding supernova debris 
are decelerated by a reverse shock (e.g., Chevalier 1982). 
This non-radiative shock is seen as very broad, high- 
velocity Lya and Ha emission features in Space Tele- 
scope Imaging Spectrograph (STIS) data (Michael et al. 
2003, 1998; Sonneborn et al. 1998). This emission, which 
was predicted by Borkowski et al. (1997), results from 
the collisional excitation of neutral H atoms from the su- 
pernova debris crossing the shock front. Using a series 
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of long-slit STIS spectra, Michael et al. (1998) mapped 
the geometry of the reverse shock, finding it to reside 
within roughly ±30° of the equator. New ground-based 
observations reported here, supplemented by additional 
observations with STIS (Sonneborn et al. 1998; Heng et 
al., in preparation), show that the broad Lya and Ha re- 
verse shock emission has increased by a factor ^4 since 
1997. Meanwhile, the X-ray emission from SNR1987A 
has brightened at a rapidly accelerating rate (Park et al. 
2005). The blast wave interaction with the hot spots now 
dominates the X-ray emission (Zhekov et al. 2005). 

In this Letter we present new spectra from the Mag- 
ellan Observatory, which demonstrate that we can con- 
tinue to study the emission from this reverse shock with 
ground-based telescopes, despite the recent demise of 
STIS. We also discuss the possibility that the shocked 
gas will produce sufficient ionizing luminosity to pho- 
toionize hydrogen atoms in the supernova debris before 
they reach the reverse shock, and thereby suppress the 
broad Lya and Ha emission. 

2. MAGELLAN OBSERVATIONS 

We observed SN1987A during the commissioning run 
of the Low Dispersion Survey Spectrograph-3 (LDSS3) 
mounted on the Clay Telescope of the Magellan Obser- 
vatory on the evening of 2005 Feb 24 (Feb 25 UT), almost 
exactly 18 years after the supernova was first discovered. 
LDSS3 has an STA0500A 4064x4064 CCD detector. At 
red wavelengths it has a pixel scale of 0"189x 1.124 A and 
an effective spectral resolution of ~5 A (R=1300) with 
a 0'.'8 (~4.5 pixel) slit width. During the observations 
the weather was clear and the seeing was ^0'.'8. The slit 
aperture was centered on the supernova and oriented at 
P.A.=-10°, aligned with the minor axis of the equato- 
rial ring. Wavelengths and velocities are measured with 
respect to the rest wavelengths of the narrow nebular 
emission along the slit. Airglow lines and emission from 
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the surrounding H n region were subtracted by fitting 
the emission from the background sky along the slit. 

The region of the long-slit spectrum around Ha re- 
sulting from a total exposure time of 900s is shown in 
Figure la. In order to flux-calibrate the LDSS3 spec- 
trum in February 2005, we interpolated between the flux 
measured in 2004 December and 2005 May in HST/ACS 
images. We used images in the F658N filter, which in- 
cludes narrow Ha and [N n] emission from the circum- 
stellar ring and hot spots. We extracted the flux over 
the same spatial region centered on the northern half of 
the ring in both the images and the LDSS spectrum, 
and summed the flux over the F658N filter bandpass 
(AA ~50 A) in the spectrum. Tracings of the resulting 
flux-calibrated spectrum are shown in Figure lb, where 
we display the total flux within ±1", as well as separate 
extractions for the northern (blucshifted) and southern 
(redshifted) sides of the nebula. 

In order to isolate the very broad emission from the 
reverse shock, we interpolated across the narrow nebular 
emission lines from the ring and hot spots (Figure lc). 
We also removed the emission hump due to rapidly fading 
Ha emission from the central remnant by subtracting a 
Gaussian profile (had we included this emission, the total 
flux would have been 6.5% higher). Interpolating across 
He I A6680 and the red [S n] doublet was somewhat sub- 
jective; we relied on what appeared to be reverse shock 
emission at ^6700 A between the two sets of lines. This 
resulted in a red bump in the reverse shock line profile. 
This is probably the best representation of the true re- 
verse shock emission, and is used for the estimate below. 
A very conservative estimate made by interpolating lin- 
early from 6660 to 6760 A (the dotted or blue line in 
Figure lc) would have resulted in a total flux about 6% 
less than the value we quote below. 

We measure a total Ha flux for the broad reverse 
shock component of 1.37(±0.15)xl0 -13 ergs s _1 cm~ 2 
on 2005 February 25. This corresponds to the black 
tracing in Figure lc, integrated from 6200 to 6860 A, 
and is continuum-subtracted using the continuum level 
shown by the straight dashed line in Figure lc. This es- 
timate includes an adjustment of +15% for those parts 
of the reverse shock on the east and west edges of the 
ring that are excluded from the aperture. Finally, adopt- 
ing E(B - 10=0-16 (Fitzpatrick & Walborn 1990) and 
R=3.1, we multiply the observed reverse shock flux by 
a correction factor of 1.453, to estimate a dereddened 
broad Ha flux of 1.99(±0.22) x 10~ 13 ergs s" 1 cm" 2 . 

3. MASS FLUX ACROSS THE REVERSE SHOCK 

The total dereddened flux we measure at Earth corre- 
sponds to a total luminosity in the broad Ha line of about 
5.63xl0 34 ergs s _1 at a distance of ~50 kpc, or roughly 
15 Lq. From this we can infer the flux of neutral hydro- 
gen atoms across the reverse shock. For each neutral H 
atom crossing the reverse shock, roughly 1 Lya and 0.21 
Ha photons will be emitted (Michael et al. 2003). Thus, 
dividing the intrinsic Ha luminosity by the energy per 
Ha photon and multiplying by a factor of 5 gives a total 
luminosity of hydrogen atoms across the reverse shock 
of Nh — 8.9 xlO 46 s _1 , or a total hydrogen mass flux of 
roughly 2.3xl0" 3 M Q yr" 1 . Michael et al. (2003) found 
that the main emitting surface area of the reverse shock 
was within ±30° of the equator, just inside the nebular 



ring. To first order, the density of neutral H atoms in 
the debris prior to crossing the reverse shock is then 

N H t 

nH ~ 47ri? 3 s tan30° 

where Njj =8.9x 10 46 s _1 is the number of hydrogen 
atoms crossing the reverse shock, i?Rs is the present ra- 
dius of the reverse shock, and t —18 yr is the time since 
outburst. Here we have assumed that the supernova de- 
bris are in free expansion, so that Rrs = Vn a t. Tak- 
ing the observed velocity Vn a — 10 4 km s _1 , we find 
Rjis ~0.16 pc. This is roughly 80% of the radius of 
the forward shock, taken to be the radius of the ring. 
This ratio is close to the theoretically-expected value for 
self-similar expansion (Chevalier 1982). The Ha flux we 
measure implies Njj ~60-70 cm~ 3 . 

Figure 2 shows the history of Nh as inferred from STIS 
observations of broad Ha and Lya since 1997 (Sonneborn 
et al. 1998; Michael et al. 1998; Heng et al., in prep.) 
and the present observation. The estimated error bars of 
the values derived from the STIS observations are large 
mainly because the corrections for narrow emission lines 
were substantially greater in the lower dispersion STIS 
spectra. From Figure 2, we see that Nh has increased 
by a factor ^4 since the first observation of broad Lya 
from the reverse shock in 1997 (Sonneborn ct al. 1998). 

4. WILL PREIONIZATION SHUT OFF THE REVERSE 
SHOCK EMISSION? 

The hot shocked gas lying immediately outside the re- 
verse shock surface is a luminous source of ionizing pho- 
tons, roughly half of which will propagate inward to pho- 
toionizc hydrogen atoms in the supernova debris before 
they reach the reverse shock. If the luminosity of these 
ionizing photons exceeds that of the hydrogen atoms, the 
broad Lya and Ha emission will vanish. Will this event 
take place, and if so, when? 

We can estimate the luminosity of ionizing photons 
from the forward shock of SNR1987A from Chandra X- 
ray observations. Spectral analysis of the LETG obser- 
vations (Zhekov et al. 2005, also Zhekov et al. in prepa- 
ration) shows that a two-shock model with temperatures 
of 0.51 and 2.7 keV, respectively, gives an excellent fit 
to the X-ray spectrum at t = 17.5 years. For such a 
model, most of the ionizing photons are EUV photons 
having energies well below the 0.4 - 10 keV Chandra 
band. From the two-shock model spectrum, we estimate 
an inward luminosity of ionizing (> 13.6 eV) photons of 
Fi ~ 3.7 x 10 45 s _1 , where we have included a factor of 
1.4 to account for the brightening of X-rays that took 
place between t = 17.5 and t = 18.0 years, and a factor 
0.5 to account for the fact that only half of the ioniz- 
ing photons will propagate inward. It follows, then, that 
the ratio of ionizing photon luminosity to hydrogen atom 
luminosity at t = 18.0 years is R t — Fi/N H ^ 0.04. 

We regard this estimate as conservative because we 
know that the complex system of shocks in SNR 1987A 
must have velocities ranging from ~ 150 km s^ 1 (as ob- 
served in the line profiles of the optical hot spots; Pun 
et al. 2002) to - 300 - 1700 km s^ 1 (as observed in the 
X-ray line profiles; Zhekov et al. 2005). We estimated Fj 
based on the best model fit to the observed X-ray spec- 
trum which requires shocks of ~ 600 and ~ 1400 km 
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s , respectively. This may be a significant underes- 
timate, because shocks having velocities in the range 
150 < Vs < 500 km s _1 might contribute substantially 
to Fi but relatively little to the Chandra band. 

At present, the X-ray luminosity of SNR1987A is in- 
creasing by a factor of ~1.7 every year (Park et al. 2005), 
more rapidly than the luminosity of hydrogen atoms 
across the shock. If present trends continue, as indi- 
cated in Figure [3 F{ should overtake Ni by about 2012 
to 2014. This event could occur 1-2 years later if the 
expanding SN ejecta have a high He abundance, because 
absorption by He atoms could reduce the effective value 
of Fi. We regard the prediction in Figure 2 as very con- 
servative, however. As we have suggested, we may have 
underestimated the value of Fi at 18 years. Moreover, 
Luo et al. (1994) predicted that, once the blast wave 
envelops the circumstellar ring, the ionizing luminosity 
should rise very rapidly to a value Fi ~ 2 x 10 48 s _1 . 
This rapid rise, and the consequent vanishing of the re- 
verse shock emission, could happen anytime during the 
next several years. 

Additional uncertainty arises from the extrapolation of 
Nh- There are two hydrodynamic scenarios that might 
bracket its evolution. The first is the self-similar solution 
describing the expansion of a supernova envelope with 
density law p{v, t) oc t~ 3 v~ n into a uniform circumstellar 
medium (Chevalier 1982). In this scenario, the radii of 
the blast wave and the reverse shock both increase as 

R x £(n-3)/n and the 

mass luminosity increases as Nh oc 
t. We regard this behavior as a probable lower limit 
to the rate of increase of Nh- But Nh must begin to 
increase more rapidly in the near future. The blast wave 
is now overtaking the dense hotspots on the circumstellar 
ring, and each such encounter will send a reflected shock 
inward toward the reverse shock. The reflected shocks 
will eventually merge with the reverse shock and bring 
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it nearly to a halt. Thereafter, Nh oc t^ n ~ 4 \ or oc t 5 
for a typical value n = 9 (Eastman & Kirshner 1989). 
We regard this as a probable upper limit to the future 
evolution of Nh- These two limiting behaviors of Nh are 
shown as the dashed curves in Figure 

5. DISCUSSION 

When the suppression of the broad Ha emission is seen, 
it will provide vital insights into the radiation hydro- 
dynamics of the developing SNR. For example, if the 
suppression occurs early, as we might expect for reasons 
noted above, it would provide evidence for substantial 
ionizing radiation from shocks too slow to contribute sub- 
stantially to the observed X-ray spectrum. Also, we don't 
expect the suppression to take place uniformly around 
the reverse shock. X-ray images (Park et al. 2005) cer- 
tainly show a substantial amount of structure. Even 
if the supernova debris have cylindrical symmetry, we 
would expect the broad Ha and Lya emission to vanish 
first near the brighter sources of ionizing radiation (i.e. 
the bright X-ray knots). With ground-based telescopes, 
such variation should be evident in the evolution of the 
line profiles. 

Yogi Berra said, "It's tough to make predictions, espe- 
cially about the future." That statement certainly ap- 
plies to the future evolution of SNR 1987A. Even so, we 
think that the reverse shock emission is likely to vanish 
soon enough that it merits continued vigilance with both 
ground- and space-based telescopes. Yogi Berra also said, 
"You can observe a lot just by watching." 
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Fig. 1. — LDSS3 spectra of SN1987A showing the Ho reverse shock emission, (a) Long-slit spectrum showing narrow emission lines 
from the ring and hot spots as well as very broad (±10 4 km s — x ) Ha emission from the reverse shock, (b) Extracted spectra. The top 
tracing (thick black line) is the total emission integrated over ±1" in the top panel, whereas the thin tracing (blue in the online edition) 
is the blueshifted side of the reverse shock at 0—1", and the dotted tracing (red in the online edition) is the redshifted side of the shock at 
-1-0". (c) Emission from the reverse shock with the narrow emission components removed. The thick black tracing shows the spectrum 
used to estimate the pure reverse shock emission. The dashed tracing (red in the online edition) includes the Ho emission from the central 
remnant, which would increase the estimated flux by +6.5%. The straight line segment (blue in the online edition) shows a conservative 
approximation for the reverse shock flux interpolated over He I A6680 and the [S II] doublet, which would be 6% less than the total flux in 
the solid tracing. The dashed black line shows the continuum level we chose in measuring the total reverse shock flux. 
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Fig. 2. — The luminosity of H-atoms vs. ionizing photons at the reverse shock in SNR 1987A. The X-ray fluxes used to derive the ionizing 
photon curve are taken from Park et al. (2005) and then converted into photon luminosity using the same two-shock model as explained 
in the text. The dotted line represent the increase of a factor 1.7 yr _1 , extrapolating from the behavior of the X-ray data over the past 
few years. The Ha data from this paper (solid dot) and from Heng et al. (in prep.; diamonds), and Lye? data from Sonneborn et al. (1998; 
unfilled circle) are used to calculate the H-atom luminosity (see text). The dashed and dash-dotted lines mark the likely lower limit (oc t) 
and upper limit (oc t 5 ) to the future evolution of the H-atom luminosity. 



Draft version February 3, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



THE REVERSE SHOCK OF SNR1987A AT 18 YEARS AFTER OUTBURST 1 

Nathan Smith 2,3 , Svetozar A. Zhekov 4 ' 8 , Kevin Heng 4 , Richard McCray 4 , Jon A. Morse 5 ' 6 , and Mike Gladders 7 

Draft version February 3, 2008 

ABSTRACT 

We use low-dispersion spectra obtained at the Magellan Observatory to study the broad Ha emission 
from the reverse shock of the infant supernova remnant SNR1987A. These spectra demonstrate that 
the spatio-kinematic structure of the reverse shock can be distinguished from that of the circumstellar 
ring and hotspots, even at ground-based spatial resolution. We measure a total dereddened Ha flux 
of 1.99(±0.22)xl0 -13 ergs s -1 cm~ 2 at an epoch 18.00 years after outburst. At 50 kpc, the total 
reverse shock luminosity in Ha is roughly 15 Lq, which implies a total flux of neutral hydrogen atoms 
across the reverse shock of 8.9xl0 46 s _1 , or roughly 2.3xl0 -3 Mq yr _1 . This represents an increase 
by a factor ~ 4 since 1997. Lyman continuum radiation from gas shocked by the forward blast wave 
can ionize neutral hydrogen atoms in the supernova debris before they reach the reverse shock. If the 
inward flux of ionizing photons exceeds the flux of hydrogen atoms approaching the reverse shock, 
this pre-ionization will shut off the broad Lya and Ha emission. The observed X-ray emission of 
SNR1987A implies that the ratio of ionizing flux to hydrogen atom flux across the reverse shock is 
presently at least 0.04. The X-ray emission is increasing much faster than the flux of atoms, and if 
these trends continue, we estimate that the broad Lya and Ha emission will vanish in <7 years. 
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1. INTRODUCTION 

The collision between the ejecta of SN1987A and its 
circumstellar ring is now in full bloom, signaling the birth 
of the supernova remnant SNR1987A. This interaction 
was predicted (Luo & McCray 1991; Luo et al. 1994; 
Chevalier & Dwarkadas 1995; Borkowski et al. 1997) 
shortly after the discovery of the circumstellar ring. It 
began about a decade after the outburst with the dis- 
covery of the first of many "hot spots" (Sonneborn et al. 
1998; Michael et al. 2000; Pun et al. 2002). These hot 
spots are thought to result when the forward blast wave 
encounters and transmits radiative shocks into dense pro- 
trusions or "fingers" pointing to the interior of the ring. 
Since then, hot spots have encircled the entire ring (Sug- 
erman et al. 2002). 

Behind the blast wave, the expanding supernova debris 
are decelerated by a reverse shock (e.g., Chevalier 1982). 
This non-radiative shock is seen as very broad, high- 
velocity Lya and Ha emission features in Space Tele- 
scope Imaging Spectrograph (STIS) data (Michael et al. 
2003, 1998; Sonneborn et al. 1998). This emission, which 
was predicted by Borkowski et al. (1997), results from 
the collisional excitation of neutral H atoms from the su- 
pernova debris crossing the shock front. Using a series 
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of long-slit STIS spectra, Michael et al. (1998) mapped 
the geometry of the reverse shock, finding it to reside 
within roughly ±30° of the equator. New ground-based 
observations reported here, supplemented by additional 
observations with STIS (Sonneborn et al. 1998; Heng et 
al., in preparation), show that the broad Lya and Ha re- 
verse shock emission has increased by a factor ^4 since 
1997. Meanwhile, the X-ray emission from SNR1987A 
has brightened at a rapidly accelerating rate (Park et al. 
2005). The blast wave interaction with the hot spots now 
dominates the X-ray emission (Zhekov et al. 2005). 

In this Letter we present new spectra from the Mag- 
ellan Observatory, which demonstrate that we can con- 
tinue to study the emission from this reverse shock with 
ground-based telescopes, despite the recent demise of 
STIS. We also discuss the possibility that the shocked 
gas will produce sufficient ionizing luminosity to pho- 
toionize hydrogen atoms in the supernova debris before 
they reach the reverse shock, and thereby suppress the 
broad Lya and Ha emission. 

2. MAGELLAN OBSERVATIONS 

We observed SN1987A during the commissioning run 
of the Low Dispersion Survey Spectrograph-3 (LDSS3) 
mounted on the Clay Telescope of the Magellan Obser- 
vatory on the evening of 2005 Feb 24 (Feb 25 UT), almost 
exactly 18 years after the supernova was first discovered. 
LDSS3 has an STA0500A 4064x4064 CCD detector. At 
red wavelengths it has a pixel scale of 0"189x 1.124 A and 
an effective spectral resolution of ~5 A (R=1300) with 
a 0'.'8 (~4.5 pixel) slit width. During the observations 
the weather was clear and the seeing was ^0'.'8. The slit 
aperture was centered on the supernova and oriented at 
P.A.=-10°, aligned with the minor axis of the equato- 
rial ring. Wavelengths and velocities are measured with 
respect to the rest wavelengths of the narrow nebular 
emission along the slit. Airglow lines and emission from 
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the surrounding H n region were subtracted by fitting 
the emission from the background sky along the slit. 

The region of the long-slit spectrum around Ha re- 
sulting from a total exposure time of 900s is shown in 
Figure la. In order to flux-calibrate the LDSS3 spec- 
trum in February 2005, we interpolated between the flux 
measured in 2004 December and 2005 May in HST/ACS 
images. We used images in the F658N filter, which in- 
cludes narrow Ha and [N n] emission from the circum- 
stellar ring and hot spots. We extracted the flux over 
the same spatial region centered on the northern half of 
the ring in both the images and the LDSS spectrum, 
and summed the flux over the F658N filter bandpass 
(AA ~50 A) in the spectrum. Tracings of the resulting 
flux-calibrated spectrum are shown in Figure lb, where 
we display the total flux within ±1", as well as separate 
extractions for the northern (blucshifted) and southern 
(redshifted) sides of the nebula. 

In order to isolate the very broad emission from the 
reverse shock, we interpolated across the narrow nebular 
emission lines from the ring and hot spots (Figure lc). 
We also removed the emission hump due to rapidly fading 
Ha emission from the central remnant by subtracting a 
Gaussian profile (had we included this emission, the total 
flux would have been 6.5% higher). Interpolating across 
He I A6680 and the red [S n] doublet was somewhat sub- 
jective; we relied on what appeared to be reverse shock 
emission at ^6700 A between the two sets of lines. This 
resulted in a red bump in the reverse shock line profile. 
This is probably the best representation of the true re- 
verse shock emission, and is used for the estimate below. 
A very conservative estimate made by interpolating lin- 
early from 6660 to 6760 A (the dotted or blue line in 
Figure lc) would have resulted in a total flux about 6% 
less than the value we quote below. 

We measure a total Ha flux for the broad reverse 
shock component of 1.37(±0.15)xl0 -13 ergs s _1 cm~ 2 
on 2005 February 25. This corresponds to the black 
tracing in Figure lc, integrated from 6200 to 6860 A, 
and is continuum-subtracted using the continuum level 
shown by the straight dashed line in Figure lc. This es- 
timate includes an adjustment of +15% for those parts 
of the reverse shock on the east and west edges of the 
ring that are excluded from the aperture. Finally, adopt- 
ing E(B - 10=0-16 (Fitzpatrick & Walborn 1990) and 
R=3.1, we multiply the observed reverse shock flux by 
a correction factor of 1.453, to estimate a dereddened 
broad Ha flux of 1.99(±0.22) x 10~ 13 ergs s" 1 cm" 2 . 

3. MASS FLUX ACROSS THE REVERSE SHOCK 

The total dereddened flux we measure at Earth corre- 
sponds to a total luminosity in the broad Ha line of about 
5.63xl0 34 ergs s _1 at a distance of ~50 kpc, or roughly 
15 Lq. From this we can infer the flux of neutral hydro- 
gen atoms across the reverse shock. For each neutral H 
atom crossing the reverse shock, roughly 1 Lya and 0.21 
Ha photons will be emitted (Michael et al. 2003). Thus, 
dividing the intrinsic Ha luminosity by the energy per 
Ha photon and multiplying by a factor of 5 gives a total 
luminosity of hydrogen atoms across the reverse shock 
of Nh — 8.9 xlO 46 s _1 , or a total hydrogen mass flux of 
roughly 2.3xl0" 3 M Q yr" 1 . Michael et al. (2003) found 
that the main emitting surface area of the reverse shock 
was within ±30° of the equator, just inside the nebular 



ring. To first order, the density of neutral H atoms in 
the debris prior to crossing the reverse shock is then 

N H t 

nH ~ 47ri? 3 s tan30° 

where Njj =8.9x 10 46 s _1 is the number of hydrogen 
atoms crossing the reverse shock, i?Rs is the present ra- 
dius of the reverse shock, and t —18 yr is the time since 
outburst. Here we have assumed that the supernova de- 
bris are in free expansion, so that Rrs = Vn a t. Tak- 
ing the observed velocity Vn a — 10 4 km s _1 , we find 
Rjis ~0.16 pc. This is roughly 80% of the radius of 
the forward shock, taken to be the radius of the ring. 
This ratio is close to the theoretically-expected value for 
self-similar expansion (Chevalier 1982). The Ha flux we 
measure implies Njj ~60-70 cm~ 3 . 

Figure 2 shows the history of Nh as inferred from STIS 
observations of broad Ha and Lya since 1997 (Sonneborn 
et al. 1998; Michael et al. 1998; Heng et al., in prep.) 
and the present observation. The estimated error bars of 
the values derived from the STIS observations are large 
mainly because the corrections for narrow emission lines 
were substantially greater in the lower dispersion STIS 
spectra. From Figure 2, we see that Nh has increased 
by a factor ^4 since the first observation of broad Lya 
from the reverse shock in 1997 (Sonneborn ct al. 1998). 

4. WILL PREIONIZATION SHUT OFF THE REVERSE 
SHOCK EMISSION? 

The hot shocked gas lying immediately outside the re- 
verse shock surface is a luminous source of ionizing pho- 
tons, roughly half of which will propagate inward to pho- 
toionizc hydrogen atoms in the supernova debris before 
they reach the reverse shock. If the luminosity of these 
ionizing photons exceeds that of the hydrogen atoms, the 
broad Lya and Ha emission will vanish. Will this event 
take place, and if so, when? 

We can estimate the luminosity of ionizing photons 
from the forward shock of SNR1987A from Chandra X- 
ray observations. Spectral analysis of the LETG obser- 
vations (Zhekov et al. 2005, also Zhekov et al. in prepa- 
ration) shows that a two-shock model with temperatures 
of 0.51 and 2.7 keV, respectively, gives an excellent fit 
to the X-ray spectrum at t = 17.5 years. For such a 
model, most of the ionizing photons are EUV photons 
having energies well below the 0.4 - 10 keV Chandra 
band. From the two-shock model spectrum, we estimate 
an inward luminosity of ionizing (> 13.6 eV) photons of 
Fi ~ 3.7 x 10 45 s _1 , where we have included a factor of 
1.4 to account for the brightening of X-rays that took 
place between t = 17.5 and t = 18.0 years, and a factor 
0.5 to account for the fact that only half of the ioniz- 
ing photons will propagate inward. It follows, then, that 
the ratio of ionizing photon luminosity to hydrogen atom 
luminosity at t = 18.0 years is R t — Fi/N H ^ 0.04. 

We regard this estimate as conservative because we 
know that the complex system of shocks in SNR 1987A 
must have velocities ranging from ~ 150 km s^ 1 (as ob- 
served in the line profiles of the optical hot spots; Pun 
et al. 2002) to - 300 - 1700 km s^ 1 (as observed in the 
X-ray line profiles; Zhekov et al. 2005). We estimated Fj 
based on the best model fit to the observed X-ray spec- 
trum which requires shocks of ~ 600 and ~ 1400 km 
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s , respectively. This may be a significant underes- 
timate, because shocks having velocities in the range 
150 < Vs < 500 km s _1 might contribute substantially 
to Fi but relatively little to the Chandra band. 

At present, the X-ray luminosity of SNR1987A is in- 
creasing by a factor of ~1.7 every year (Park et al. 2005), 
more rapidly than the luminosity of hydrogen atoms 
across the shock. If present trends continue, as indi- 
cated in Figure [3 F{ should overtake Ni by about 2012 
to 2014. This event could occur 1-2 years later if the 
expanding SN ejecta have a high He abundance, because 
absorption by He atoms could reduce the effective value 
of Fi. We regard the prediction in Figure 2 as very con- 
servative, however. As we have suggested, we may have 
underestimated the value of Fi at 18 years. Moreover, 
Luo et al. (1994) predicted that, once the blast wave 
envelops the circumstellar ring, the ionizing luminosity 
should rise very rapidly to a value Fi ~ 2 x 10 48 s _1 . 
This rapid rise, and the consequent vanishing of the re- 
verse shock emission, could happen anytime during the 
next several years. 

Additional uncertainty arises from the extrapolation of 
Nh- There are two hydrodynamic scenarios that might 
bracket its evolution. The first is the self-similar solution 
describing the expansion of a supernova envelope with 
density law p{v, t) oc t~ 3 v~ n into a uniform circumstellar 
medium (Chevalier 1982). In this scenario, the radii of 
the blast wave and the reverse shock both increase as 

R x £(n-3)/n and the 

mass luminosity increases as Nh oc 
t. We regard this behavior as a probable lower limit 
to the rate of increase of Nh- But Nh must begin to 
increase more rapidly in the near future. The blast wave 
is now overtaking the dense hotspots on the circumstellar 
ring, and each such encounter will send a reflected shock 
inward toward the reverse shock. The reflected shocks 
will eventually merge with the reverse shock and bring 
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it nearly to a halt. Thereafter, Nh oc t^ n ~ 4 \ or oc t 5 
for a typical value n = 9 (Eastman & Kirshner 1989). 
We regard this as a probable upper limit to the future 
evolution of Nh- These two limiting behaviors of Nh are 
shown as the dashed curves in Figure 

5. DISCUSSION 

When the suppression of the broad Ha emission is seen, 
it will provide vital insights into the radiation hydro- 
dynamics of the developing SNR. For example, if the 
suppression occurs early, as we might expect for reasons 
noted above, it would provide evidence for substantial 
ionizing radiation from shocks too slow to contribute sub- 
stantially to the observed X-ray spectrum. Also, we don't 
expect the suppression to take place uniformly around 
the reverse shock. X-ray images (Park et al. 2005) cer- 
tainly show a substantial amount of structure. Even 
if the supernova debris have cylindrical symmetry, we 
would expect the broad Ha and Lya emission to vanish 
first near the brighter sources of ionizing radiation (i.e. 
the bright X-ray knots). With ground-based telescopes, 
such variation should be evident in the evolution of the 
line profiles. 

Yogi Berra said, "It's tough to make predictions, espe- 
cially about the future." That statement certainly ap- 
plies to the future evolution of SNR 1987A. Even so, we 
think that the reverse shock emission is likely to vanish 
soon enough that it merits continued vigilance with both 
ground- and space-based telescopes. Yogi Berra also said, 
"You can observe a lot just by watching." 
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Fig. 1. — LDSS3 spectra of SN1987A showing the Ho reverse shock emission, (a) Long-slit spectrum showing narrow emission lines 
from the ring and hot spots as well as very broad (±10 4 km s — x ) Ha emission from the reverse shock, (b) Extracted spectra. The top 
tracing (thick black line) is the total emission integrated over ±1" in the top panel, whereas the thin tracing (blue in the online edition) 
is the blueshifted side of the reverse shock at 0—1", and the dotted tracing (red in the online edition) is the redshifted side of the shock at 
-1-0". (c) Emission from the reverse shock with the narrow emission components removed. The thick black tracing shows the spectrum 
used to estimate the pure reverse shock emission. The dashed tracing (red in the online edition) includes the Ho emission from the central 
remnant, which would increase the estimated flux by +6.5%. The straight line segment (blue in the online edition) shows a conservative 
approximation for the reverse shock flux interpolated over He I A6680 and the [S II] doublet, which would be 6% less than the total flux in 
the solid tracing. The dashed black line shows the continuum level we chose in measuring the total reverse shock flux. 
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Fig. 2. — The luminosity of H-atoms vs. ionizing photons at the reverse shock in SNR 1987A. The X-ray fluxes used to derive the ionizing 
photon curve are taken from Park et al. (2005) and then converted into photon luminosity using the same two-shock model as explained 
in the text. The dotted line represent the increase of a factor 1.7 yr _1 , extrapolating from the behavior of the X-ray data over the past 
few years. The Ha data from this paper (solid dot) and from Heng et al. (in prep.; diamonds), and Lye? data from Sonneborn et al. (1998; 
unfilled circle) are used to calculate the H-atom luminosity (see text). The dashed and dash-dotted lines mark the likely lower limit (oc t) 
and upper limit (oc t 5 ) to the future evolution of the H-atom luminosity. 



